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Abstract: To date, the role of dopamine Dl receptor (DRDl) polymoiphism in schizophrenia 
remains controversial. We earned out a meta-analysis to determine whether DiiD7 polymorphism 
influences the risk of schizophrenia. We examined whether rs4532 and rs5326 genetic variants 
are related to the etiology of schizoplirenia, using a meta-analysis. Relevant case-control studies 
were retrieved by database searching and selected according to established inclusion criteria. 
A total of ten studies were identified and included in our meta-analysis, nine for rs4532, with 
1,941 cases and 2,480 controls, and four for rs5326, with 1,285 cases and 1,195 controls. No 
significant association was found between the rs4532 locus and schizophrenia. For the rs5326 
locus, the guanine-adenine (GA) genotype was associated with schizophrenia as a risk factor (for 
GA vs guanine-guanine [GG], odds ratio [OR] =1.36, 95% confidence intei'val [CI]: 1.15-1.61, 
P<0.001). The GA genotype of rs5326 increased the risk of schizophrenia, but there was no 
association between rs4532 and schizophrenia. These data may provide references for case- 
control studies in schizophrenia in future. 
Keywords: susceptibility, SNP, genetics, psychiatry 

Introduction 

Schizophrenia is a common mental disorder affected by the interaction of muhiple 
genetic and environmental factors, with a median incidence of 15.2/100,000 persons 
worldwide.' Evidence from family, adoption, and twins studies supports that its heri- 
tability is up to 80%.^ Many efforts have been made over the years, in order to clarify 
the pathological mechanisms, but a consistent conclusion hasn't been achieved given 
that nongenetic factors, such as social and psychological factors, have an impact on the 
occurrence of schizophrenia and increase the difficulty of the identification of suscep- 
tible genes. ^ Nevertheless, genome- wide association studies have identified a genetic 
association by comparing common single nucleotide polymorphisms across the human 
genome,^ and a number of copy number variants have been suggested as susceptibility 
factors for schizophrenia.* To date, dysregulation of dopaminergic neurotransmission 
has been implicated in the pathogenesis of schizophrenia,^ '' and the reduction of pre- 
frontal cortical dopamine neurotransmission likely leads to schizophrenia.' 

Dopamine Dl receptor {DRDl) gene polymorphisms roughly play a role in the occur- 
rence of schizophrenia by affecting the expression of the Dl receptor. Linkage associa- 
tion between DRDl and schizophrenia has been reported in Chinese,^ American,' and 
Portuguese populations.'" Increasing evidence indicates that the polymorphism in the 5' 
untranslated region of DRDl is associated with the mediation of the therapeutic response 
to antipsychotics," Moreover, present studies'^"''' show that deficits in the DRDl transmis- 
sion are associated with the cognitive and negative symptoms of schizophrenia and that 
there is a significant decrease in DRDl expression in the basal ganglia of schizophrenic 
patients. '"^ In contrast, some studies have reported that there was no genetic association 
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between the DRDl gene and the risk of schizophrenia."'"'* 
Therefore, the role of the DRDl gene in schizophrenia 
remains controversial, and a consistent conclusion has yet 
to be achieved. The human DRDl gene rs4532 (A-48G) and 
rs5326 (G-94A) loci are two commonly studied polymorphisms 
that are located on the 5' UTR region. The Rs4532 locus has 
been associated with therapeutic response to antipsychotics," 
psychotic symptoms, and prefrontal cortex-related cognitive 
performance.'-^" As for the rs5326 locus, this may act as a 
short-open-reading frame preceding the receptor cistron and 
encode a small peptide that inhibits receptor translation.^' Thus, 
it is necessary and sufficient to estimate the role of rs4532 and 
rs5326 loci in the pathogenesis of schizophrenia. 

There are several possible explanations for this conflicting 
evidence, such as small sample size, ethnic differences, and 
publication bias. To systematically synthesize these scattered 
studies, a potentially useful method is to use meta-analytic 
techniques, which can combine the results of several studies 
to generate a single estimate of the major effect with enhanced 
precision.^^ Results from multiple studies are incorporated in 
an unbiased fashion.^^ Meta-analysis is becoming important 
in psychiatric genetics because of its great improvement and 
rapid increase of the number and size of data sets. Conse- 
quently, in our present study, we carried out a meta-analysis 
of studies examining the association between DRDl gene 
polymorphism and the risk of schizophrenia. 

Methods 

Identification and eligibility 
of relevant studies 

To identify studies eligible for our meta-analysis, three online 
electronic databases (PubMed, Web of Science, and China 
National Knowledge Infrastructure [CNKI]) were searched 
(the last search update was December 2013). The follow- 
ing key words were used in the literature search: DRDl, 
dopamine receptor 1, dopamine Dl receptor, dopamine 
receptor Dl, schizophrenia, polymorphism, and variation. 
The reference lists from identified articles and potentially 
relevant review articles were also screened to identify addi- 
tional relevant studies. Studies included in this meta-analysis 
met the following inclusion criteria: (1) used a case-control 
design; (2) included patients with schizophrenia and healthy 
controls; (3) offered available allele or genotype frequencies; 
(4) the most recent of studies with the same or overlapping 
data published by the same authors; and (5) published in 
an English language, peer-reviewed journal. Major reasons 
for exclusion of studies were: (1) no control population; 
(2) duplicate of earlier publication; and (3) no usable geno- 
type frequency data. 
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Data extraction 

According to the inclusion criteria listed above, two of the 
authors extracted information fi'om all eligible publications 
independently. If there were missing data, we contacted the 
responsible authors to retrieve additional data not contained 
in the original paper. Any disagreement was resolved through 
discussion until the two authors reached a consensus. The fol- 
lowing data were included from each study: the first author's 
last name, publication year, region, race, study design, source 
of controls, diagnostic criteria, detection method, and num- 
bers of genotypes between cases and controls. 

Statistical analysis 

The strength of the association between two loci and schizo- 
phrenia was measured by odd ratios (OR), with 95% confi- 
dence intervals (CI). Pooled effect sizes across studies were 
performed by a random effects model, which evaluates the 
likely effect size across different populations and takes het- 
erogeneity across studies into account. The random effects 
model is different from a fixed effects model, which evaluates 
the most likely effect size from multiple studies, hypothesiz- 
ing that they are sampled fi^om a single population, but it can 
be biased by high heterogeneity across studies. 

The degree of heterogeneity between studies was per- 
formed using a g-statistic.^^-^'' A P-value >0.05 for the g-test 
indicated a lack of heterogeneity, and P- value <0.05 indi- 
cated an existence of heterogeneity. The P is the proportion 
of observed variance in effect sizes attributable to the true 
differences among studies. A conventional interpretation of F 
is that it defines bounds for low (<25%), moderate (-50%), 
and high (>75%) heterogeneity. ^"'^^ Subgroup analysis was 
performed by type of region. 

An estimate of publication bias was performed by visual 
inspection of a funnel plot, in which the standard error of 
log (OR) of each study was plotted against its log (OR). An 
asymmetric plot indicated a possible publication bias, and the 
degree of asymmetry was tested using Egger's test (/'<0.05 
was considered to indicate significant publication bias).^* 

Sensitivity analysis was carried out to assess the poten- 
tial influences of any single study on the pooled effect size. 
This was performed by sequentially omitting each study 
within each meta-analysis, in order to check for significant 
alterations to the pooled effect sizes. All statistical tests were 
two-sided. The meta-analysis was performed using Stata ver- 
sion 10.0 (StataCorp LP, College Station, TX, USA). 

Results 

The search strategy retrieved 44 potentially relevant studies. 
According to the inclusion criteria, ten articles were identified 
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Table I Baseline characteristics of the studies included in meta-analysis 



Author 


Year 


Region 


Race 


Study design 


Source of controls 


Diagnostic criteria 


Detection method 


Cichon et aP' 


1994 


German 


Caucasian 


Case-control 


Population-based 


DSM-III-R 


PCR-RFLP 


Dmitrzak et aP° 


2006 


Poland 


Caucasian 


Case-control 


Population-based 


DSM-IV/ICD-IO 


PCR-RFLP 


Dollfus et al" 


1996 


France 


Caucasian 


Case-control 


Population-based 


DSM-III-R 


PCR-RFLP 


Du et al" 


201 1 


People's Republic 
of China 


Mongolian 


Case-control 


Population-based 


DSM-IV 


Taqman®* 


Hu et al" 


2005 


People's Republic 
of China 


Mongolian 


Case-control 


Population-based 


CCMD-III 


PCR-RFLP 


Iwata et al™ 


2003 


Japan 


Mongolian 


Case-control 


Population-based 


DSM-IV 


PCR-RFLP 


Kojima et aP' 


1999 


Japan 


Mongolian 


Case-control 


Hospital-based 


DSM-IV 


PCR-RFLP 


Zhang et al" 


2010 


People's Republic 
of China 


Mongolian 


Case-control 


Hospital-based 


DSM-IV 


Taqman®* 


Zhu et al" 


201 1 


People's Republic 
of China 


Mongolian 


Case-control 


Hospital-based 


DSM-IV 


MALDI-TOF 


Zhang et al" 


2010 


People's Republic 
of China 


Mongolian 


Case-control 


Hospital-based 


DSM-IV 


Taqman®* 



Note: *Taqman SNP genotyping assay. 

Abbreviations: CCMD, Chinese Classification and Diagnostic Criteria of Mental Disorders; DSM, American Diagnostic and Statistical Manual of Mental Disorders: ICD. 
International Classification of Disease; MALDI-TOF, matrix-assisted laser desorption/lionization time of flight mass spectrometry; PCR, polymerase chain reaction; RFLP, 
restriction fragment length polymorphism; SNP, single nucleotide polymorphism. 



and included in our meta-analysis (Table 1). Among the ten 
eligible studies, nine studies, with 1,941 cases and 
2,480 controls, were included in the meta-analysis for the 
rs4532 locus and four studies, with 1,285 cases and 
1,195 controls, in the analysis for rs5326 locus in Table 2. 

Meta-analysis result 

Table 3 lists the main results of pooled OR for the loci of 
rs4532 and rs5326 and schizophrenia. No significant asso- 
ciation was found between the rs4532 locus and schizophre- 
nia (for adenine-guanine [AG] vs adenine-adenine [AA] 
[OR =0.76, 95% CI: 0.54-1. 07, P=0.1 15] or guanine-guanine 
[GG] vs AA [OR =0.90, 95% CI: 0.64-1.26,^=0.526]). For 
the rs4532 locus in the allele contrast model (guanine [G] 
vs adenine [A]), OR =0.95, 95% CI: 0.76-1.20, /'=0.665. 
For the rs5326 locus, the GA genotype was associated with 



schizophrenia, as a risk factor (for GA vs GG, OR =1.36, 
95% CI: 1 . 1 5-1 .6 1 , P<0.00 1 ), while no association between 
the AA genotype and schizophrenia was observed (for AA 
vs GG, OR =0.89, 95% CI: 0.5 1-1 .54, P=0.675). For rs5326 
locus in the allele contrast model (A vs G), OR = 1 .24, 95% 
CI: 1.16-1.32, P<OMl. The A allele was associated with 
schizophrenia as a risk factor. The forest plot is shown in 
Figure 1. 

Subgroup analysis 

Examination of heterogeneity among the rs4532 locus 
group across studies was significant. However, the genotype 
frequency of rs4532 locus had diverse distribution among 
different races. We carried out the subgroup analysis, by 
race, in the rs4532 locus group. Races were divided into a 
Mongolian subgroup (including Chinese and Japanese) and 



Table 2 Genotype frequencies of rs4532 and rs5326 of all studies included in the meta-analysis 



First author 


Year 


rs4532 (case/control) 




HWE test 
in control 




rs5326 (case/control) 




HWE test 
in control 


AA 


AG 


GG 




P-value 


GG 


GA 


AA 




P-value 


Cichon et al" 


1994 


4/5 


16/20 


16/20 


0 


1.000 












Dmitrzak et aP° 


2006 


64/51 


186/489 


157/159 


1 3 1 .272 


<0.00l 












Dollfus et al" 


1996 


1/4 


14/40 


47/1 17 


0.069 


0.793 












Du et aP" 


201 1 












85/58 


72/32 


5/1 1 


3.692 


0.055 


Hu et aP' 


2005 


90/134 


21/52 


6/2 


1.559 


0.212 












Iwata et aP° 


2003 


39/126 


8/18 


I/I 


0.161 


0.688 












Kojima et aP' 


1999 


1 14/126 


33/22 


I/O 


0.954 


0.329 












Zhang et al" 


2010 


280/265 


88/108 


5/6 


1.819 


0.177 


209/23 1 


144/124 


20/24 


1.725 


0.189 


Zhu et al" 


201 1 


299/252 


81/92 


5/6 


0.532 


0.466 


204/216 


150/1 14 


31/20 


0.909 


0.341 


Zhang et al" 


2010 


275/256 


85/103 


5/6 


1.455 


0.228 


18/24 


I4I/I 15 


206/226 


3.066 


0.080 



Note: All P<0.05 were considered to be statistically significant. 
Abbreviations: A, adenine; G, guanine; HWE, Hardy-Wei nberg equilibrium. 
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Table 3 Main 


results of pooled ORs for rs4532 and rs5326 










OR 


95% CI 


P-value 


P' 




VJ 1 IILJIIInLIVJII IIICLIIUU 


rs4532 A/G 












Random effect model 


AG vs AA 


0.76 


0.54-1.07 


0.1 15 


71.9% 


<0.00l 


Random effect model 


GG vs AA 


0.90 


0.64-1.26 


0.526 


0.0% 


0.633 


Random effect model 


Mo+oo vs r\f\ 


0.81 


0.62-1.05 


0.1 13 


57.8% 


0.0 1 5 


IXdllUUIII CllCl_L IIIUUCl 


rsS326 G/A 












Random effect model 


GA vs GG 


1.36 


I.I5-I.6I 


<0.00l 


0.0% 


0.947 


Random effect model 


AA vs GG 


0.89 


0.5 1- 1.54 


0.675 


59.2% 


0.061 


Random effect model 


GA+AA vs GG 


1.29 


1. 10-1.52 


0.002 


0.0% 


0.8879 


Random effect model 



Notes: ^/^ represents the variation in OR attributable to heterogeneity. ''P-value of Q-test for the heterogeneity test. 
Abbreviations: CI, confidence interval; OR, odds ratio; A, adenine; G, guanine; h, heterogeneity. 



Caucasian subgroup (including German, Polish, and French). 
Table 4 lists the results of the subgroup analysis. No signifi- 
cant association was found between the rs4532 locus and 
schizophrenia in the Mongolian and Caucasian subgroups. 
Examination of heterogeneity across studies was also not 
significant, which showed that the heterogeneity among the 
rs4532 locus group was generated by different races. 

Sensitivity analysis 

A sensitivity analysis was carried out for each meta-analysis 
to assess the influence of every single study. There were no 



significant differences in the corresponding pooled ORs when 
removing one study at the time from each meta-analysis, 
indicating that our results were stable and reliable. 

Publication bias 

A funnel plot was performed to assess publication bias of 
the literature. The funnel plot is shown in Figure 2. Egger's 
test was used to supply the statistical evidence for funnel 
plot symmetry. For rs4532 locus, the results did not show 
any evidence of publication bias (for AG vs AA [t=0.82, 
P=0.439] or for GG vs AA [1=235, P=0.051]). For the rs5326 



study 
ID 



rs4532 AG vs AA 



% 

OR (95% CI) Weight 



rs4532 GG vs AA 



% 

OR (95% CI) Weight 



Cichon et aF — 
Dmitrzak et aP" — 

Dollfus et al'^ 

Du et aP 
Hu et aP 
Iwata et al™ 
Kojima et aP' 
Ztiang et al" 
Ztiu et aP3 
Ztiang et al^^ 

Overall (f =71 .9%, P=0.000) 



1.00(0.23. 4.35) 
0.30 (0.20, 0.45) 
-> 1.40(0.14, 13.61) 
0.60(0.34, 1.07) 
1.44(0.58, 3.56) 
1.66 (0.91, 3.01) 
0.77(0.56, 1.07) 
0.74(0.53, 1.05) 
0.77(0.55, 1.07) 
0.76 (0.54, 1.07) 



4.12 
14.71 

I. 98 
12.14 
7.96 

II. 80 
15.86 
15.65 
15.77 
100.00 



1.00 (0.23, 4.35) 
0.79 (0.51, 1.21) 
1.61 (0.17, 14.75) 
4.47 (0.88, 22.63) 
- 3.23 (0.20, 52.86) 
-) 3.31 (0.13, 82.18) 
0.79 (0.24, 2.61 ) 
0.70(0.21,2.33) 
0.78 (0.23, 2.57) 
0.90 (0.64, 1.26) 



5.26 

61.81 

2.31 

4.32 

1.46 

1.10 

7.91 

7.91 

7.91 

100.00 



Study 
ID 



rs5326 GA vs GG 



% 

OR (95% CI) Weight 



rs5326 AA vs GG 



OR (95% CI) Weight 



Du et a\^ 
Zhang et al" 
Zhu et aP3 
Zhang et aP^ 

Overall (P=0.0%, P=0.947) 



o 



-) 1.54(0.90, 2.62) 10.00 ^ 

1.28(0.95,1.74) 30.73 

1.39(1.02,1.90) 2969 

1.35(0.99,1.83) 2958 

1.36(1.15,1.61) 100.00 



0.382 



2.62 



Figure I Forest plot for association between rs4532 and rs5326 loci and schizophrenia. 
Note: Weights are from random effects analysis. 

Abbreviations: CI, confidence interval; OR, odds ratio; A, adenine; G, guanine. 



0.31(0.10,0.94) 15.76 

0.92(0.49.1.72) 27.96 

1.64(0.91,2.97) 28.87 

0.82(0.43,1.56) 27.40 

0.89(0.51,1.54) 100.00 
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Table 4 Subgroup analysis 


between rs4532 and the risk of schizophrenia 








Race 


Surf 


msry of pooled ORs 




Heterogeneity test 




OR 


9S% CI 


P-va.1 u e 


12a 
f 


rfi 


Mongolian 












AG vs AA 


0.84 


0.66-1.07 


0.157 


41.5% 


0.128 


GG vs AA 


I.IO 


0.60-2.04 


0.761 


1.3% 


0.408 


A<j+tjtj VS AA 


0.87 


0.69-1.10 


U.Z JO 


4 1 Q°/ 
*i 1 . 7/0 


u. 1 zo 


Caucasian 












AG vs AA 


0.52 


0.19-1.42 


0.205 


48.4% 


0.144 


GGvs AA 


0.82 


0.55-1.23 


0.338 


0.0% 


0.795 


AG+GG vs AA 


0.55 


0.29-1.04 


0.066 


21.4% 


0.280 



Notes: ^/^ represents the variation in OR attributable to heterogeneity. ^P-value of Q-test for the heterogeneity test. 
Abbreviations: CI, confidence interval: OR, odds ratio; A, adenine: G, guanine: h, heterogeneity. 



locus, the results did not show any evidence of publication 
bias (for GA vs GG [^=1.94, P=0.192] of for GA vs GG 
[^=-2.31,P=0.147]). 

Discussion 

This was a meta-analysis testing the association between the 
genetic polymorphisms of the DRDl gene loci of rs4532 and 
rs5326 and the risk of schizophrenia. In our study, the GA 
genotype of the rs5326 locus increased the risk of schizophre- 
nia, and the A allele was also associated with schizophrenia 



as a risk factor, while no association was observed between 
rs4532 locus and schizophrenia. 

Dopamine receptors play a pivotal role in regulating emo- 
tion, cognition, and motor behavior, as well as in neuroendo- 
crine signaling. Prefrontal dopamine D 1 receptor signaling 
is required for temporal control and for guiding intentional 
movements.'^ Abnormal modulation of dopamine receptor 
signaling and disorder of dopaminergic nerve function may 
cause a series of neuropsychiatric defects.^* The DRDl gene 
has been considered a potential susceptibility gene for various 



rs4532 AG vs AA 



rs4532 GG vs AA 



O 
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O 
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O 

G) 
O 

4- 

o 

S 
L_ 

<s> 

■a 
k. 
m 
■o 
c 
m 

(O 



0.2 0.4 

log OR 

Figure 2 Funnel plot for evaluating the publication bias in schizophrenia. 
Abbreviations: A, adenine; G, guanine; OR, odds ratio. 
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neuropsychiatric illnesses.' "-^* Furthermore, several previous 
studies have reported that the DRDl gene was associated 
with alcohol dependence,^' hyperactivity disorders,"" nicotine 
addiction,'" ''^ and other behavioral disorders. 

The DRDl gene rs4532 and rs5326 loci are located in 
the 5' untranslated region, where the nucleotide sequence is 
involved in transcription but is not being translated into amino 
acids. A previous meta-analysis'*^ indicated that the DRDl 
gene might contain a genuine susceptibility allele, for which 
the ORforrs4532 reached 1.18 (95% CI: 1.01-1.38). Further- 
more, rs4532 was associated with antipsychotic treatment in 
schizophrenia," " and also to working memory and cognitive 
behavior related to the prefrontal cortex.^" These results sug- 
gest that rs4532 likely has a potential link with the etiology of 
schizophrenia. Rs4532, due to a G to A substitution, may act 
as a short open reading frame preceding the receptor cistron 
and encode a small peptide that inhibits receptor translation.^' 
However, the positive or negative association between the 
rs5326 locus and schizophrenia hasn't been determined.'^-''' 

Heterogeneity is usually a concern in a meta-analysis. 
In our study, evidence of heterogeneity was observed in the 
rs4532 group. This might have been partially due to differ- 
ent races because the genotype frequencies changed greatly 
in different regions. When subgroup analysis was carried 
out to explain the heterogeneity, we found that differences 
in races could lead to the heterogeneity among the studies. 
The funnel plot was not completely symmetrical, which 
indicated the potential presence of publication bias. This 
could have resulted from the following causes: firstly, nega- 
tive results were less often published compared with positive 
results, which made it hard for us to collect them for analysis; 
second, the literature collection wasn't comprehensive, and 
the number of published studies in our meta-analysis wasn't 
sufficiently large; and third, studies written in other languages 
were difficult to collect and analyze. 

Several limitations in our study should be noted. First, 
the sample sizes of several studies included were fairly small 
and lacked a strong statistical power. Furthermore, the fac- 
tors of sex'* and age may affect the incidence of the disease. 
In addition, we didn't consider the role of different environ- 
ments, and it is possible that environment may also influence 
the risk of schizophrenia. Further, epigenetic mechanisms 
altering the chromatin structure, such as histone acetyla- 
tion and deoxyribonucleic acid methylation, may mediate 
effects of environmental factors in transcriptional regula- 
tion of specific genes and may be a prominent factor in 
gene-environmental interaction.'''' Thus, further studies are 
needed to investigate the effect of genetic susceptibility and 



environment in different ethnicities and between males and 
females, to fully attest to our results. In the future, research 
into the role of DRDl gene polymorphisms in schizophrenia 
endophenotypes could be useful. 

In conclusion, our meta-analysis demonstrated that the 
A allele of the DRDl rs5326 polymorphism might be a poten- 
tial risk factor in schizophrenia. As the number of eligible 
studies was limited in this meta-analysis, these results need 
further investigation. 
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